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Abstract. The static and high-frequency differential conductivity of a one- 
dimensional superlattice with parabolic miniband, in which the dispersion law is 
assumed to be parabolic up to the Brillouin zone edge, are investigated theoretically. 
Unlike the earlier published works, devoted to this problem, the novel formula for the 
static current density contains temperature dependence, which leads to the current 
maximum shift to the low field side with increasing temperature. 

The high-frequency differential conductivity response properties including the 
temperature dependence is examined and opportunities of creating a terahertz 
oscillator on Bloch electron oscillations in such superlattices are discussed. 

Analysis shows that superlattices with parabolic miniband dispersion law may be 
used for generation and amplification of terahertz fields only at very low temperatures 

(r-.o). 



PACS numbers: 72.10, 72.60 
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1. Introduction 

In present work, we study theoretically the static and high-frequency conductivity of a 
semiconductor superlattice (SL). Unlike the earlier published numerous works, devoted 
to this problem, where the conventional cosine-type model was used for the conduction 
miniband, here a dispersion law is considered in form of a truncated parabola, i. e. the 
dispersion law is assumed to be parabolic up to the Brillouin zone edge. Such a problem 
statement is of interest, among others, from the view point of opportunities of creating 
a terahertz oscillator on Bloch electron oscillations in SLs. In works different 
variants of realization of such opportunity were discussed and it was mentioned that 
the main obstacle consists in using the non-optimal SL structures, in particular, the SL 
with cosine-type miniband. 

Thus the theoretical investigations of electric properties of SLs with other dispersion 
laws are necessary, all the more so since the modern technology allows to vary widely 
the form of the potential relief and the SL energy spectrum. 

The main condition for realization of Bloch oscillator consists in existence of 
negative high-frequency differential conductivity on that regions of current-voltage 
characteristic where the static differential conductivity is positive. In [2] it was shown 
that this condition holds, in particular, in SL with parabolic miniband. But this result 
was obtained in the limiting case T — >■ 0. Here we find the temperature dependence of 
conductivity of such SL and define the temperature criterion by which the mentioned 
condition holds practically. 

This article is structured as follows. In Section [2] we derive an expression for 
static conductivity of SL with parabolic miniband, which is valid for any temperatures. 
In Section [3] we derive the corresponding expression for high-frequency differential 
conductivity. Section H] presents the conclusions of our work. 

2. Static distribution function and current-voltage characteristic 

The electron energy in the SL lowest miniband is [1] 



where p is quasimomentum, d is SL period, x axis being directed along the SL axis, 
£{p±) is in-plane electron energy, n'^h'^ /md"^ = A is double miniband width, m is effective 
electron mass. 

In quasi-classical situation (A ^ eEd, h/r, where r is electron momentum 
relaxation time, e is electron charge), the current density in electric field E*°*(t) may be 
found by solving Boltzmann equation with collision integral within r-approximation: 



where -fo(p) is equilibrium electron distribution function, F{p,t) is unknown 
distribution function perturbed due the electric field. Below we use dimensionless 
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variables by changing pd/{iTh) -> p, E^°^edT/{nh) E*°*, T/A T, t/r t (T 
is temperature in energy units). 

With the field E*°*(t) is directed along the SL axis (E*°*(t) = {E^°\t), 0, 0)), we have 
^(P,^) = fo{p±)fiP,t), ^o(p) = /o(p±)/o(p), where /o(p) is equilibrium distribution 
function, normalized to the carrier density n (/o(p±) being normalized to unity). Thus, 
the function f{p,t) satisfies the following equation 



fo{p)-f{p,t), {-l<p<l) 



dt dp 

with periodicity conditions /(l,t) = /(— l,t). 

In a static field E^°^{t) = E = const, and denoting f{p) = fc{p, E, T), we get 

E^ = fo-fe, {-l<p<l). 
We consider non-degenerate electron gas, so that 



(3) 



(4) 



fo{p,T) = 2n 



-1 



exp - 



p^ 
2T 



y2^erf ( 

where erf(z) is error function. In the low temperature limit (T 
reduces to the function used in [Ij: go{p) = 2n6{p). 

The exact solution of (jlj) with periodicity condition, /c(— - 
form [5] 



(5) 

0) the relation ([5]) 

= /c(l), takes the 
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1 < p <1. 



(6) 



In limiting case E ^ reduces to (JSj). In another limiting case, T — > 0, we get 
the distribution function found in [1]: 



\([l-expi-2/E)]-\ 0<p<l, 



g{p, ) = ^ exp ^ [g^p^2/E) -1< p < 



(7) 



The function ([6]) satisfies the same normalization condition as the equilibrium 
function /o 

-j Up,E,T)dp = n (8) 
-1 

and, therefore, it makes the integral of right-hand side of formula (jl]) vanish. Besides, 
the integral of left-hand side of the Boltzmann equation (jlj) vanishes too, because of the 
periodicity condition mentioned. The distribution function fc{p,E,T) at several values 
of E and T is shown in figure [H 
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Figure 1. Distribution function fc{p) at various values of the driving field and 
temperature. 1) E ^ 1, T ^ 0.1; 2) E ^ 0.1, T = 0.1; 3) £; = 2, T = 2; 4) 
E = 1, T = 0.005. The dashed curve 5 represents function g{p) at i? = 1. 



The current density j in the direction of SL axis can be found (in dimensional units) 
by a conventional way 

Trh/d 



J 



ed 



2'7ihm 



pfcip) dp. 



(9) 



-Trh/d 



By substitution function ([6]) into ([9]) we get 

-1 



j{E,T) = E 



2erf 



1 



X 



erf 



sinh 

'2T J 
( VT 1 



exp 
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\2E^ 
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erf 



\E^/2 
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(10) 



Here j is expressed in units of jo = nel^dj-nh,^ while all the quantities are written in 
dimensionless form. 

Equation (fTO!) determines the current- volt age characteristic for the parabolic 
miniband SL with the current density temperature dependence taking into account. 
To warrant numerical stability we present formula ffTOj) in the following form 



2{E,T) = Ea{E,T), a{E,T) = l- 
where a{E,T) is the conductivity and 

E' 



I 2 exp(-0.5/T) + A(E,T) 



erf(l/V2T) 



l/E 



A{E,T) 



Tsinh(l/E) 



s^E^' 
2T 



exp 



s sinh s ds. 



(11) 
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The value of A{E, T) can be estimated numerically with high accuracy. 
Expanding the exponent in a power series we get 

ME,T).lt^^^, (13) 
where functions Gn{E) are defined by recurrent formula 

Go = E coth (^-^) - E^, Gn = Go + 2nE^ [{2n + - 1] . (14) 

As Gn{E) G [0, 1/ (2n + 3)), series ( fT3l) converges quickly. As numerical experiments 
show, first four terms of series ( |T3l) give good approximation at T > 0.5. 

At \E\ — ^ we have A{E, T) 0, so in low fields {\E\ <^ 1) in linear approximation 
on E we have 



where angle brackets mean averaging over the equilibrium distribution. Note that 
the conductivity temperature dependence in low fields (the expression within round 
brackets in f[T^ ) is close to the analogous dependence for the miniband cosine model 
(/i(l/2T)//o(l/2T), In{z) being the modified Bessel function). 
In high fields {\E\ > 1) we have 

^ ' E^TVTeTiil/V2T)^o (2^)" ' " h W\{2n - 2k + 5)V ^ ' 
where are Bernoulli numbers. 

For low temperatures (T <^ 1), using f[T^ . we get 



a(E,T) 



erf (1/V2T 



i5sinh(l/E)^"P 12^;+!/^^"^-^ 



(17) 



As numerical experiments show, formula flTTI) gives good approximation at T < 0.07. 
In limiting case T ^ from ( |T7I) we get the expression that was found in [1] 

j=j{E) = E- . ] (18) 
smh(l/_fc) 

From (fT5][T6|) it follows that j ~ E at |E| < 1 and j ~ 1/E at |E| > 1. Therefore 
at fixed temperature T = fix the function j{E,T) reaches its maximum at some value 
E = EciT) > and negative differential conductivity is realized aX E > Ec{T) (see 
figure [2]). 

Note, that EciT) decreases with increasing temperature. Essentially, that Ec value 
does not depend on the temperature at all in the cosine model: Ec = 1/11^ 0.318. 

The parametric representation of dependence EciT) is defined by equation ad = 0, 
where ad = dj/dE is the differential conductivity. Using fllllll2p . we get 

cr,(E, T) = 1 + ( \e coth f - t] [aiE, T) - 1] - 

(19) 
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Figure 2. Current-voltage characteristic at different values of temperature. 1) T = 0; 
2) T = 0.01; 3) T = 0.1; 4) T = 0.5; 5) T = 1. 



E 



Thus function Ec{T) is defined implicitly by equation 
1 



^•'^^erf 



TA{E,T) = Ecoth 



T) + exp 



1 

2T 



(20) 



2 

and it is sufficient to solve this equation at E > 0. 

The numerical solution of equation ( l20l) at versus T is presented in figure [3l 
Note that dependence E = Ec(T) is monotone so the inverse function Tq = Tc{E) 

exists. To investigate behavior of function Tc{E), consider first the case of high 

temperatures T ^ 1. Expanding all functions in a power series on 1/T and neglecting 

all terms o(l/T^), we get 

(45E^ + 22.5^2 + 1.8) tanh^ (1 /E) - {36E + 3E) tanh (1 /E) - 9E^ - 1.5 



c 



By that 



(9^2 + 4)tanr {1/E) - 6^ tanh (1/^) - 3 



lim TciE) = +CX), 



where E\ ^ 0.29104955 is the root of equation 

(9E^ + 4) tanh^ f— ) - 6Etanh [— 



EJ ^-^~--V£- ^ °- 
Consider now the case of low temperatures T <^ 1. Using ([1 



we get 



Tc{E) ^ 2E^ 
By that 



E'^ tanh ( — ^ sinh ^— 



lim TciE) = 0, 



- 1 



1-2E tanh 



E 



(21) 
(22) 
(23) 

(24) 
(25) 
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Figure 3. The dependence E = Ec{T). The dashed curve Ec = I/tt represents Ec 
for cosine model. 



where E2 ~ 0.373681745 is the root of equation 
.2 - / 1 \ . - / 1 



Enanh y—j sinh j = 1. (26) 

Therefore function Ec{T) is defined for T > and 

hm EciT) = E2, hm Ec(T) = Ei. (27) 

T^O T~»+oo 

3. High-frequency differential conductivity 

In this section we will determine the induced superlattice current in the presence of an 
external electric field given by 

E^°\t) = E + Eq cos ut, (28) 

where u is measured in units of r~^. Within the scope of quasi-classical conditions the 
value of E is arbitrary. Assuming the amplitude of variable field Eq to be much smaller 
then the static field E, consider the time-dependent field in linear approximation. The 
distribution function may be found in a form 

f{p, E, T, t) = f,{p, E, T) + /i(p, E, T, iu) expi-tcot), (29) 

here fi{p,E,u!) satisfies the following equation [2] 

E^ + il- = -^0^ (30) 
with periodicity condition fi{—l,E,uj) = fi{l,E,uj) and by 

fj\{p,E,T,u;)dp = 0. (31) 
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It is easily to show that required solution is 

00 ' E 



fi{p,E,T,uj) 



Up,E,T) + fJp 



E 



T 



(32) 



1 — iuj, 

With the help of (132|) the dynamic (high-frequency) differential conductivity can be 
found by a conventional way. The result is 



ai{E,T,uj) = — 
uoE 



j{E,T)-j 



E 



.T 



(33) 



.1 — zcu 

From fl55]) it follows that at — the value ai{E,T,uj) tends to static differential 
conductivity (fT9l) 

lima^{E,T,u;) = a^{E,T). (34) 
Using (fTOj) . we get 



Re ai{E, T, u) 



2ujE 



J 



E 



1 + iu 



.T 



J 



E 



1 — iuj 



T 



(35) 



For numerical computations we present expression (135!) in a form 

1 



Re o"i(i?, T, u) 



1 
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sinh2(l/E) + sin2(cj/E)' 


-1 


ttT cuEerf (l/y2T) 
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exp 
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At T 



from fl36|) we get the expression presented in [2] 

1 cosh(l/E)sin(a;/i?) 

ite (Ti(£/, U, 00) = 



oo^ 



ooE 



sinh2(l/E) + sm^ioo/E) 



(36) 



(37) 



The opportunities of creating a terahertz oscillator on Bloch electron oscillations 
in SLs are defined by conditions of existence of negative high-frequency differential 
conductivity on that regions of current- volt age characteristic where the static differential 
conductivity is positive [2]-|l]. These conditions would prevent development of 
undesirable domain instabilities (Gunn effect). 

Let Q = eEd/h be the Bloch oscillations frequency which in normalized 
measurement units is equal to irE. Then the static differential conductivity is positive 
at n < and negative at > fic, where fic = ttEc{T) e (0.914, 1.174). Thus the 
conditions of low-frequency domain instability suppression are defined by that values of 
parameters oo and Q, for which 

ai{Q,T,oo) < 

The existence of such conditions for regarded model of dispersion law was discovered 
in [2] in limiting case T — > 0. 



(3^ 
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a) b) 

Figure 4. The regions of negative high-frequency differential conductivity at 
parameter plane {n,uj). a.) T ^ 0, ^ 1.174. b) T = 0.01, f^c = 1.06. 



But conditions fl55]) prove to be very sensitive to temperature increasing. In figure H] 
the regions in parameter space {Q, T, u) are presented in which the high-frequency 
differential conductivity is negative. 

The boundary hues of these regions are defined by condition Reai{E,T,uj) = 0. 
At these hues we have uj ~ kQ, k = 1,2,.... Thus the frequencies at which the high- 
frequency differential conductivity changes sign are multiples of the Bloch frequency. 

Note the existence of regions of low-frequency domain instability suppression at 
T = and absence of such regions at T = 0.01. 

The dependence of function cti on parameters Q and uj is presented in figure [51 

Note that by temperature increasing the oscillations of di become suppressed at 
Q < Qc and negative high-frequency differential conductivity disappears. 

4. Conclusion 

In present paper, an exact distribution function has been found of the carriers in 
the lowest parabolic miniband of a SL, placed in the dc electric field, parallel to SL 
axis. The novel formula for the static current density in SL contains temperature 
dependence, which leads to the current maximum shift to the low field side with 
increasing temperature. 

We have obtained explicit expression for high-frequency differential conductivity 
at arbitrary temperature. It was shown that high-frequency differential conductivity is 
very sensitive to temperature of SL. We have compared high-frequency electron behavior 
at different temperatures and exhibited the drastic change in the character of regions 
where the high-frequency differential conductivity is negative. In particular we have 
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a) 



b) 



Figure 5. a) Driving field dependence of tiigh-frequency differential conductivity at 
u; = 10. 1) T = 0, = 1-174; 2) T = 0.01, = 1-163; 3) T = 0.1, = 1-06. 
b) Dependence of high-frequency differential conductivity onajatfi^l. 1)T = 0; 
2) T — 0.01; 3) T — 0.1. At such temperatures the static differential conductivity 
'^d = cTiL^o is positive. 



discovered that the possibihty of low-frequency domain instabihty suppression may be 
reahzed only at T — * 0. 

In summary, our analysis shows that SLs with parabolic miniband dispersion law 
may be used for generation and amplification of terahertz fields only at very low 
temperatures (T < 0.01 A). 

The numerical estimations of the effects predicted are reduced, in general, to 
measurement units of electric field and temperature. At d = 10^^ cm, r = 10~^^ s, 
A ^ 10^^ eV we get that units for E and T are ~ 2 ■ 10^ and ~ 100 K respectively. 
Thus the condition T < O.OIA is equivalent to T < 1 K. 
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